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Density Functional Theory (DFT) is one of the most widely used methods for "ab initio" calculations 
of the structure of atoms, molecules, crystals, surfaces, and their interactions. Unfortunately, the 
customary introduction to DFT is often considered too lengthy to be included in various curricula. 
An alternative introduction to DFT is presented here, drawing on ideas which are well-known from 
thermodynamics, especially the idea of switching between different independent variables. The central 
theme of DFT, i.e. the notion that it is possible and beneficial to replace the dependence on the external 
potential v{r) by a dependence on the density distribution n(r), is presented as a straightforward 
generalization of the familiar Legendre transform from the chemical potential jj, to the number of 
particles A^. This approach is used here to introduce the Hohenberg-Kohn energy functional and to 
obtain the corresponding theorems, using classical nonuniform fluids as simple examples. The energy 
functional for electronic systems is considered next, and the Kohn-Sham equations are derived. The 
exchange-correlation part of this functional is discussed, including both the local density approximation 
to it, and its formally exact expression in terms of the exchange-correlation hole. A very brief survey 
of various applications and extensions is included. 



I. INTRODUCTION 



The predominant theoretical picture of solid-state 
and/or molecular systems involves the inhomogeneous 
electron gas: a set of interacting point electrons mov- 
ing quantum-mechanically in the potential field of a 
set of atomic nuclei, which are considered to be static 
(the Born-Oppenheimer approximation). Solution of 
such models generally requires the use of approximation 
schemes, of which the most basic — the independent 
electron approximation, the Hartree theory and Hartree- 
Fock theory — are routinely taught to undergraduates 
in Physics and Chemistry courses. However, there is an- 
other approach — Density Functional Theory (DFT) — 
which over the last thirty years or so has become in- 
creasingly the method of choice for the solution of such 
problems (see Fig. |l|.). This method has the double ad- 
vantage of being able to treat many problems to a suf- 
ficiently high accuracy, as well as being computationally 
simple (simpler even than the Hartree scheme). Despite 
these advantages it is absent from most undergraduate 
and many graduate curricula with which we are familiar. 

We believe that this omission stems in part from the 
tendency of the existing books and review papers on 
DFT, e.g. Refs. to follow the historical path of 

development of the theory. Although appropriate for 
a thorough treatment, this approach unnecessarily pro- 
longs the introduction and grapples with problems which 
are not directly relevant to the practitioner. It is our pur- 
pose here to give a brief and self-contained introduction 
to density functional theory, assuming only a first course 
in quantum mechanics and in thermostatistics. We break 
with the traditional approach by relying on the analogy 
with thermodynamics B. In this formulation, the use 
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FIG. 1. One indicator of the increasing use of DFT is the 
number of records retrieved from the INSPEC databases by 
searching for the keywords "density", "functional" and "the- 
ory" . This is compared here with a similar search for keywords 
"Hartree" and "Fock", which parallels the overall growth of 
the INSPEC databases (for any given year, approximately 
0.3% of the records have the Hartree-Fock keywords). 



of the density distribution as a free variable arises in a 
natural manner, as do more advanced concepts which are 
central to recent developments in the theory ||^, e.g. the 
exchange-correlation hole and generalized compressibili- 
ties. The discussion is sufficiently detailed to provide a 
useful overview for the beginning practitioner, and the 
relatively novel point of view may also prove illuminat- 
ing for those experienced researchers who are not familiar 
with it. We hope that the availability of such an introduc- 
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tion will encourage teachers to include a one or two hour 
class on DFT in courses on quantum mechanics, atomic 
and molecular physics, condensed matter physics, and 
materials science. 

The general theoretical framework of DFT, involving 
the Hohenberg-Kohn free energy FiiK[n{r)], is presented 
in Sec. II, which for simplicity focuses on classical sys- 
tems. The generalization to the quantum-mechanical 
electron gas is given in Sec. Ill, together with the dis- 
cussion of the Kohn-Sham equations and of the local 
density approximation, which is the simplest practical 
approximation for the exchange-correlation energy. Var- 
ious issues relating to the accuracy of this approach are 
discussed in Sec. IV, followed by a summary in Sec. V. 

II. GENERAL THEORY 

In this section, a unified treatment of thermodynamics 
and density functional theory is presented. For simplic- 
ity, the case of a classical interacting system of point par- 
ticles will be discussed first. Although classical DFT has 
its own applications, e.g. liquids Q, the reader is advised 
to keep in mind electronic systems, which will be the sub- 
ject of the next section. Thus, Eqs. (|l]) through (||), to be 
derived here using classical notation, are equally applica- 
ble to quantum-mechanical systems, where Hilbert space 
with its position and momentum operators replaces the 
classical phase space and its scalar coordinates. 

A. Thermodynamics: a reminder 

We begin by rederiving the equations of thermodynam- 
ics from statistical mechanics [Q. Consider a classical 
system of M interacting particles in a container of vol- 
ume V. The many-body Hamiltonian is: 

T^MB = T + U , (1) 

where T = J2iLi pf /2m is the kinetic energy, and U — 
J2i<j '^il'^^i ~ ^]\) is the interaction energy, assuming a 
simple pair potential u{r). Here and p; are the posi- 
tions and momenta of the particles, and m is their mass. 
We consider the grand-canonical ensemble, where the 
system is in contact with a heat reservoir of tempera- 
ture T and a particle reservoir with chemical potential /i. 
It is well-known from statistical physics that the grand 
potential, which is the free energy in this case, is given 
by: 

n{fi,T,v) = -riogs, (2) 

where S is the grand partition function, 

-(/^'^'^)-f:]i^Trexp(-^^M^) , (3) 




FIG. 2. (a) The Legendre transform which gives F{N) 
corresponds to describing the curve by the properties 

of its tangents: minus their slopes A'' = —(dil/dn) and their 
intercepts with the energy axis, F — n + fiN. The fact that 
the derivative AF/ AN is equal to /i follows from asking the 
question: if two neighboring lines intercept the energy axis at 
a distance AF from each other, and have slopes which differ 
by AA'', how far from the axis will they cross each other? 
(b) The Legendre transform back from F{N) to ^l{fi) has a 
similar interpretation. The minimization suggested in Eq. (^) 
corresponds to studying a family of lines with a fixed slope 
fi, which pass through points (A'', F) on the free energy curve. 
Their intercepts, — F — fiN, have a minimum (marked by 
an asterisk) for that line which is tangent to the curve. 

the temperature is in energy units (i.e. = 1), and 
the classical trace, Tr, represents the 6M-dimensional 
phase-space integral (the division by M\ compensates for 
double counting of many-body states of indistinguishable 
particles). 

It follows directly from these definitions that the ex- 
pectation value of the number of particles in the sys- 
tem is given by a derivative of the grand potential, 
TV = (M) = -{dVt/dix). The convexity of the ther- 
modynamic potential |8| implies that iV is a monoton- 
ically increasing function of /i. Other partial derivatives 
of give the values of additional physical quantities, 
such as the entropy, S = —{dVL/dT) and the pressure 
P = —{dfl/dV). This may be summarized by writing 
dn = -Ndn- SdT - PdV. 

A basic lesson of thermodynamics is that in different 
contexts it is advantageous to use different ensembles. 
For example, in studying systems where the number of 
particles rather than the chemical potential is fixed, it is 
preferable to use the Helmholtz free energy ||^, which 
is obtained from the grand potential r2 by a Legen- 
dre transform: F{N,T,V) = n{fi{N),T,V) + ^iiN)N. 
Here fJ-{N) is no longer an independent variable, but 
a function of N obtained by inverting the relationship 
N = N{^i,V,T) = -(dVt/dn). The derivative of F 
with respect to the "new" free variable N is equal to 
the "old" free variable /j,. The derivatives with respect 
to the other variables are unchanged (but are taken at 
constant N rather than at constant ji). We thus write 
dF = fidN - SdT - PdV. 

For the purpose of comparison with DFT, it is useful to 
make a variation on the inverse Legendre transform which 
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expresses in terms of F, and to define the following 
"grand potential function" , which depends explicitly on 
both fj, and N: 



n^{N,T,V) = F{N,T,V)-fiN 



(4) 



This function gives the original grand potential of Eq. (Q) 
when minimized with respect to N, i.e. when the deriva- 
tive (dF/dN) — /i vanishes, which is equivalent to the 
condition N ~ N{fi, T, V) conventionally used in the in- 
verse Legendre transform. For other values of N, the 
function n^,{N,T, V) describes a "cost" in free energy of 
having a configuration with the "wrong" number of elec- 
trons. For a geometric interpretation of Legendre trans- 
forms, including the minimization procedure of Eq. (|^), 
see Fig. H 



B. Nonuniform systems and the Hohenberg— Kohn 
theorem 



The discussion above can be generalized in a quite 
straightforward manner to the treatment of particles in 
an external potential v(r). The many-body Hamiltonian 



Hmb = T + V+U , 



(5) 



where the potential energy, V — ^i^i): been 

added. The grand potential and the partition function 
are defined as before, Eqs. (§) and (||), but they now 
depend on the potential function v{r) rather than on the 
scalar volume V. In this sense, 51 — Jl T, [ii(r)]) is 
now a functional |0 of v{r) as well as a function of fi 
and T — the square brackets denote functional variables 
(it is also implicitly a functional of the pair potential 
u{r)). As is well known, the potential v{r), is an energy 
which is measured from an arbitrary origin, i.e. shifting 
the potential by a constant does not affect the physics 
of the system. It is convenient here to set this origin at 
the chemical potential, i.e. to take ^ — 0. Equivalently, 
one may define the new functional variable as v{r) — /x, 
as Jl depends only on this difference, and not on v and fi 
separately pT[ |. 

The functional derivative of 51 with respect to the new 
variable gives the density distribution of the particles, 
n{r) = {p{r)) = Si}/5v{r), where p{r) = Y^fti '^('' ~ ^i) 
is the unaveraged density. Using a (functional) Legendre 
transform as above, we can define a new free energy which 
depends on n{r) rather than on u(r), and is called the 
Hohenberg-Kohn free energy: 



-FHK[^^(r)] = ^ / drn{r)v{r) , 



(6) 



where the explicit temperature variable has been omit- 
ted, and u(r) on the right hand side is chosen to cor- 
respond to the given n{r) (that such a choice is possi- 
ble follows from the "generalized convexity" of the free 



energy |^,|T2|). The partial and functional derivatives of 
-FHK['T-(r)] are given by the usual rules for Legendre trans- 
forms: d^HK = — 5dr— / dr v{r)Sn{r). 

The direct generalization of the free energy function of 
Eq. (0) is the free energy functional: 



fly[n{r)] = FHK[n(r)] + / drn{r)v{r) , 



(7) 



with w(r) and n(r) treated as independent functional 
variables. If this free energy functional is minimized with 
respect to n(r) at constant v{r) (and given T, etc.), the 
relation 
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HK 



Sn{r) 



(8) 



is obtained. For n(r) and v{r) obeying this physical re- 
lation, the free energy functional is equal to the grand 
potential by inspection. The existence of a functional of 
n(r) with this property is one of the basic tenets of DFT, 
and is the (second) Hohenberg-Kohn theorem jl^. 

Note that below we will use Eq. (^), which also fol- 
lows directly from the properties of Legendre transforms. 
Discussion of the Hohenberg-Kohn theorem, Eq. (^, is 
nevertheless important even in a Lcgendre-transform- 
based introduction to DFT, because the free-energy- 
minimization procedure embodied in it is central both 
to forming a physical intuitive picture of DFT, and to 
devising efficient numerical schemes for solving the DFT 
equations in practice [14[ . 



C. Examples 

We next apply the expressions above to two (simplis- 
tic) physics problems: finding the density distribution of 
air and of water in the ecosphere. 

1) Air: First we consider air, using the approximate 
model of an ideal gas. For such a gas the particles do 
not interact, u{r) = 0, and the partition function can be 
evaluated directly. The Hamiltonian, Eq. (^), reduces to, 
Ti-MB = T+V. The grand partition function, S of Eq. (||), 
may be expressed as S = X]m=o ^i*^/-^' ~ exp(Si), 
in terms of the partition sum for a single particle. Si = 
(27r?i)^'^ J drdp exp (— (p^/2m + v) /T) . The momentum 
integral is trivial, giving Si = /(dr/Aip)exp (— w(r)/r), 
where volume is normalized by the "thermal wavelength" 



At = ^jh^lmT. One finds from Eq. (|) that fl = -TSi. 
In this case the density distribution is easily found di- 
rectly: ri(r) = ^Jl/5w(r) = A^'^ exp (— u(r)/T) , but for 
pedagogical purposes we proceed to evaluate the Legen- 
dre transforms explicitly. 

Inverting the ideal gas n\v\ relationship just derived 
gives t;(r) = — Tlog(nAy), and proceeding with the func- 
tional Legendre transform of Eq. (H) gives 
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F^K[n{r)] = / drn(r)/(n(r),r) , (9) 



term. The Helmholtz free-energy per particle in the uni- 
form fluid is thus taken to be 



where 

/(n,T) = r(log(nA|)-l) , (10) 

is the free energy per particle. The latter is equal to 
the Helmholtz free energy per particle, F{N,T,V)/N, 
evaluated for a uniform ideal gas with n = N/V. 

The DFT free energy functional, Eq. (^, for an ideal 
gas is 

ny[n{r)] = J dr n(r)/ (n(r),T) + J drn(r)w(r) . (11) 

Minimizing this energy functional with respect to n{r) 
gives V = —d{fn)/dn = — Tlog(nAip), which is equiva- 
lent to n oc exp(— w/T), as discussed above. 

The density distribution of the atmosphere in the grav- 
itational field of the Earth may now be determined. The 
external potential, measured from the chemical potential, 
is v{r) = mgh — /i, where h is the height and g ~ lOm/s^ 
is the acceleration due to Earth's gravitational field (as 
we take ft, = at Earth's surface, we have subtracted /i 
explicitly here). Substituting this in the expression for 
the density, we find that the density decreases exponen- 
tially with height, n(r) cx exp(— ft//), with a length-scale 
I = T/mg. Taking m as the mass of a nitrogen molecule, 
m ~ 5 X 10~^^kg, and the temperature as T = 4 x 10~^^ 
Joules (~ 17°C), gives Z = 8 x lO'^m or 8 kilometers, less 
than the height of the Everest (in practice, the atmo- 
sphere is not in equilibrium, T ^ const., and this descrip- 
tion becomes increasingly inaccurate at higher altitudes). 

2) Water: Finding the density distribution of water re- 
quires a different analysis, because the ideal gas is a good 
description of real fluids only for low densities such as 
those found in the atmosphere. At higher densities, the 
interactions, u{r)^ must be taken into account. An ex- 
act, explicit evaluation of the partition function and its 
derivatives is no longer possible, and we are compelled to 
use models and approximations. The inhomogeneity of 
the distribution will be given an approximate treatment 
here, by assuming that the Hohenberg-Kohn energy may 
still be expressed as in Eq. (H) , in terms of the Helmholtz 
free energy per particle for a uniform system. This is 
a "local density approximation" , a concept which is in 
common use in DFT. For simplicity, we will use the van 
der Waals model for f{n,T). In this model [|l5|, both 
attractive and repulsive interactions of real atoms and 
molecules (typically included in u{r) for large and small 
r, respectively) are taken into account. Correspondingly, 
two simple modifications in the ideal-gas expressions for 
/, Eq. ( p^ are made: first, a "higher order in n" attrac- 
tive term is added, and second, an "excluded volume" 
&, representing the "hard core" of real particles, is sub- 
tracted from the volume per particle n^^ in the thermal 



f{n) = T{3 log(AT) - log(n-i - b) - l) - an , (12) 

which replaces Eq. (^0|). Correspondingly, the pressure 
is 

'' = -[dv),-^d^^-^b~''-^ (1^) 

which is a more familiar expression. 

The free-energy functional to be minimized with re- 
spect to n{r) is again that of Eq. (pT|), but the extremum 
condition, 

dinf) f b , A3 \ 

'^ = -^--^[^b^''^^b)^'''' 

(14) 

may now have several solutions, and the value of n which 
gives the lowest free energy must be selected. In the 
present model, for temperatures which are not too high 
and for a limited range of values of v, there are two com- 
peting local minima, corresponding to the liquid and the 
gas phases of the fluid ||l^ . 

We now return to water, representing a typical location 
on Earth as a unit of area (one square meter), covered 
by 10® moles of water (one mole ~ 6 x 10^^ rnolecules). 
We use the same potential v = mgh — fi and the same 
values of T and g as above. The mass of a water molecule 
is m ~ 18 X lO^^kg/niole, and we choose the values of 
a ~ 0.48Paxm^/mole^ and 6 ~ 16x 10~^m3/mole to give 
the boiling point of water as 100 degrees Celsius at a pres- 
sure of one atmosphere ~ lO^Pa, and the particle den- 
sity of the liquid water at that point as n ~ 56moles/m3, 
corresponding to a mass density of lO^kg/m^. Finding 
the density distribution, n(h), requires the inversion of 
Eq. (|lj) which is not readily available analytically; there- 
fore we have performed this exercise numerically. The 
chemical potential fi plays the role of a Lagrange multi- 
plier, imposing the constraint Jdrn(r) — N, and finding 
it requires a few trials or iterations (using the Newton- 
Raphson algorithm). 

The resulting value of ^ is —2.5385 x 10^ Joules/mole 
— several significant digits must be kept because v (or 
V — ii) varies by only 180 Joules/mole for each kilome- 
ter of altitude. The density distribution n{r) contains 
both liquid and gaseous regions. The liquid region or 
"ocean" occupies the lower altitudes, and a gaseous re- 
gion occupies higher altitudes. In the lower region, the 
fluid is more or less incompressible, with the density at 
ft = higher by only a couple of percent than that at 
"sea level" , ft ~ 1750m. Above this point one finds wa- 
ter vapor with a density of ~ 4moles/m3, corresponding 
to a pressure of lO^Pa or one-tenth of an atmosphere 
(real water deviates significantly from our van der Waals 
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model — a more accurate description of f(n) would give 
~ 2.5 X lO'^Pa). At this low density, the equation of state 
of the water vapor does not differ significantly from that 
of an ideal gas, and indeed one finds an exponential decay 
with further increase in altitude, as discussed above for 
air (the lengthscale in this case is ~ 14km, as the water 
molecules are lighter). 

Obviously, water and air coexist on the surface of 
Earth, and taking this into account would modify the 
results of our examples (most importantly, the atmo- 
sphere begins only at sea level). Such situations can be 
addressed by generalizing the free-energy density, nf, 
allowing it to depend on both densities neater and nair- 
With a sufhciently accurate parameterization of /, a good 
description of many nonuniform thermodynamic systems 
can thus be obtained. However, for such a local descrip- 
tion to hold, the spatial variations in n(r) (or equiva- 
lently, in w(r)) must be very slow on the scale of the 
range of the potential u(r). In contrast, one may be inter- 
ested in studying how the density distribution n(r) near 
a liquid-vapor interface changes gradually, on a scale of 
Angstroms, from that of the liquid to that of the gas. 
Such problems are typical applications of the DFT of 
classical systems, and must employ nonlocal model func- 
tionals 

A discussion of such nonlocal classical functionals is 
beyond the scope of the present paper, but we would like 
nevertheless to mention the following two points: (a) As 
opposed to gases, classical liquids and solids are charac- 
terized by strong and complicated correlations between 
particles, even when the interaction potential u{r) is sim- 
ple. This makes the task of finding good -FhkM func- 
tionals nontrivial. (b) The general framework of DFT 
can nevertheless give accurate descriptions of real sys- 
tems, provided that one uses good "anchoring points" . 
For example, the hard-sphere liquid has been studied 
extensively numerically, and accurate information for it 
is available. A good description of the liquid-solid tran- 
sition in several more general systems may be found by 
perturbing in the difference between the actual u(r) and 
that of the hard sphere system, Uhsif), e.g. by taking 
FhkM = Fhs + / dr(u--Uhs) {SFiiK/Su{r)) , where in the 
last term, the functional derivative is equal to 5Q/Su{r) 
which is the pair correlation function. With this as back- 
ground, we now proceed to treat electronic systems in 
the next section. 



III. APPLICATION TO ELECTRONS 

The electron is much lighter than an atom or molecule, 
and thus has a relatively large thermal wavelength, e.g. 
At — h/^/niT ~ 17A at room temperature, more than an 
order of magnitude larger than the typical inter-electron 
distance. Its treatment must be quantum-mechanical. 



with the Hamiltonian of Eq. (||) replaced by the corre- 
sponding operator, Hmb = T + V + U. The statistical- 
mechanics derivation above is affected by this, e.g. the 
trace (Tr) in the definition of the partition function, 
Eq. (||), is taken over the Hilbert space rather than 
over the classical phase space. However, the thermody- 
namic considerations and Legendre transforms, includ- 
ing Eqs. (^, (|^) and (^) for nonuniform systems, are 
unchanged. 

As electrons are Fermions, they form a degenerate gas 
at the prevailing high densities, and their energy in the 
cases of interest here does not deviate considerably from 
its ground-state value. Indeed, electronic DFT was de- 
veloped in Refs. and |17| as a ground-state theory. 
Correspondingly, we will from here on take the zero tem- 
perature limit, T 0, in all our equations and expres- 
sions (cf. Ref. |l^). In this limit the partition function 
is dominated by a single quantum-mechanical state with 
an integer number of electrons, M = N, and the grand 
potential is equal to the ground-state energy (using our 
fi — convention). The Hohenberg-Kohn free-energy, 
-FhkM, is then, by Eq. (^), equal to the total energy 
minus the potential energy, i.e. to the internal energy. 

For noninteracting quantum mechanical particles, 
u(r) = 0, the internal energy is just the kinetic energy, for 
which we introduce the notation Fni[n] = -FhkM |„(r)=o' 
Whereas in the classical case, the noninteracting or ideal- 
gas system had simple local expressions for the energy 
functionals, here the relationships between the energy, 
the density distribution, and the confinement potential 
are nontrivial and nonlocal. 

Practical implementations of DFT require an explicit 
construction of the Hohenberg-Kohn free-energy func- 
tional, -FhkM- It is customary to write ^hkM for in- 
teracting electrons as a sum of the noninteracting kinetic 
energy, Fni[n], and two interaction terms — the electro- 
static energy and the exchange-correlation energy |l9f|: 



^HK[«(r)] = i^„i[n(r)] + Ses[n(r)] + ^,c[n(r)] , (15) 

where the last term, i?xc[?^]j is defined as the remainder 
and thus contains everything that is not included in the 
first two terms. Each of the three terms on the right 
hand side is in principle a functional of the independent 
variable n(r). Only the second term — the electrostatic 
energy — is easily expressed explicitly: 



r / X, f drdr' 



n(r) n(r') , (16) 



The first and last terms are much more complicated: 
knowledge of the former implies a full understanding 
of the quantum-mechanical noninteracting problem; the 
latter contains all of the many-body physics, and is in 
principle even more complex. 

In the following, we briefly introduce two alternative 
methods for confronting this situation: (i) construction 
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of explicit approximate expressions for both -Fni['^] and 
Exc[n], and (ii) the orbital method developed by Kohn 
and Sham, which uses the noninteracting Schrodinger 
equation to evaluate Fni[n], with only the smaller term, 
Exc [n] , replaced by explicit approximations to the desired 
complicated functional. 



A. Explicit functionals 

One of the simplest models of the electronic structure 
of atoms was developed by Thomas ||2^ and Fermi ||2l| ] 
already in the late 1920's. For the nonintcracting part of 
the calculation, they used a local density approximation: 



F„i[n(r)] ~ C 



,5/3 



(r) dr , 



(17) 



with C = joi^TT ) ' - 2.87 in atomic units; here 
Cn^^^{r) is the kinetic energy density of a uniform 
electron-gas of density n(r). They also approximated 
the Coulomb interaction energy of the electrons by the 
electrostatic term only, which corresponds in the present 
language to taking -ExcM = 0. Using these simplifica- 
tions and the spherical symmetry of an atom, analytical 
progress could be made. 

Both approximations made in the Thomas-Fermi 
model are expected to be accurate at very high densities. 
The approximation made for the noninteracting part is 
accurate when the density of electrons changes slowly in 
space relative to the Fermi wavelength, or equivalently, 
if the density is sufficiently high that n{ry^^ is (much) 
larger than the spatial rate of change of the density. 
The approximate treatment of the interactions, which 
ignores exchange and correlation effects, is the leading 
order of a high-density expansion, in powers of rs/ao 
where rs = (S/inn)^^/^ is the Wigner-Seitz radius, and 
ao = h^/me^ = 0.53 A is the Bohr radius. The Thomas- 
Fermi model indeed has applications for very dense mat- 
ter, and is useful in the description of certain stars [^ . 
Unfortunately, the results for more down-to-earth sys- 
tems are rather poor. For example, it is known that in 
the Thomas-Fermi model no molecules can form — the 
dissociated atoms always have a lower energy |^ . 

The Thomas-Fermi method was extended over the 
years in two main directions. At first, approximate ex- 
pressions for the exchange-correlation energy were in- 
cluded, but this did not lead to significant improvements 
in the results. Such improvements were obtained only 
when an element of nonlocality was taken into account 
in the noninteracting kinetic energy term, i^ni[n]. This 
was achieved by including gradient terms (e.g. terms pro- 
portional to |Vnp in the integrand of Eq. (pT|), with 
n-dependent prefactors), and lead in particular to the 
possibility of modeling chemical bonds. However, the 
description of electronic structure by the Thomas-Fermi 
model and its extensions remains qualitative to date P4| . 



B. The Kohn— Sham equations 

In 1965, Kohn and Sham |]l^ made a major step to- 
wards quantitative modeling of electronic structure, by 
introducing an orbital method by which -Fni[?^] can be 
evaluated exactly. In other words, in order to evaluate 
the kinetic energy of N noninteracting particles given 
only their density distribution n(r), they simply found 
the corresponding potential, called i;efr(r), and used the 
Schrodinger equation. 



^2 \ 



(18) 



such that n(r) = X^iLi The states here are 

ordered so that the energies e.^ are non-decreasing, and 
the spin index is included in i. If cat is degenerate with 
eN+i (and also at finite temperatures l^)j fractional oc- 
cupations fi are to be used, n(r) = J2i=i /i|''/'i(r)P, but 
if only spin-degeneracy is involved, the result for the den- 
sity is not affected. The kinetic energy is then given by 
F„i[n(r)] = J2f=l{'^S^\'^^) = E^i^* - /drn(r)woff(r), 
where ti is the kinetic energy operator for the ith electron 

In practice, it is the external potential of a given sys- 
tem which is known, not the density distribution or the 
effective potential. One may find the effective potential 
by taking a functional derivative of the three-term ex- 
pression for -FHi<-[n], Eq. (|l5|), and rearranging the terms: 



(19) 



where we have used Eq. (g), SF[n\/Sn = —v, for both 
the interacting and the noninteracting system. The elec- 
trostatic potential is here 



ip{r) = 



dr 



, n(r') 



(20) 



and the exchange-correlation potential is defined as 



Wxc(r) = 



Sn{r) 



(21) 



Given a practical approximation for Exdn], one obtains 
Wxc(r), and can thus find i'cff(r) from ri(r) for a given 
v{r). 

The set of equations described above is called the 
Kohn-Sham equations of DFT, and must be solved self- 
consistently: Ucff(r) determines n{r) in Eq. (p^, and is 
determined by it in Eq. ([l9|). They provide a mecha- 
nism for minimizing the functional Ey[n\ (or r2t,[n]) of 
Eq. (^, despite the fact that i^ni['T-] is known only im- 
plicitly. This method assumes that the external poten- 
tial v{r) and the total number of electrons N are given, 
and does not use the chemical potential fi (if a nonzero 
value of n were restored, it would simply shift both sides 
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of Eq. (|T^) by a constant, and thus drop out of consid- 
eration). Together with any expHcit approximation for 
the exchange-correlation term (see below) , it provides an 
efficient scheme for finding n(r) and the ground state en- 
ergy 1 26 1 for a system of N interacting particles. Note 
that for the interacting system, only the energy and its 
derivatives (including n(r)) are accessible — the compli- 
cated many-body wavefunctions do not take part in this 
scheme. 

Historically, additional properties |^ of the Kohn- 
Sham noninteracting system, e.g. the band structure for 
crystals, have also provided surprisingly accurate predic- 
tions when compared with experiments. In fact, the 
agreement between the calculated Kohn-Sham Fermi 
surface and the measured one was so remarkable for some 
systems |Q, that it motivated analyses of soluble (per- 
turbative) models for which the difference between the in- 
teracting and noninteracting Fermi surfaces could be cal- 
culated explicitly, and shown not to vanish |^^. Clearly, 
the accuracy of DFT predictions for ground-state ener- 
gies and density distributions can be improved by find- 
ing better practical approximations for E'xcM, whereas 
improving the accuracy of such band-structure calcula- 
tions may require "going back to the drawing board" and 
devising other, more appropriate, calculational schemes 






' \ ' 
\ 

\. 




£ 


£ 






tot 


-0.1 








-0.2 




















xc 




-0.3 









1 



5 6 



r 



FIG. 3. The exchange-correlation energy per electron, 
£xc(^) (in Hartrees) of a uniform electron g&s, as a func- 
tion of the density, parameterized by the Wigner-Seitz ra- 
dius (in Bohr radii), according to the interpolation formula 
of Ref. [^. The dashed line is the corresponding total 
energy per electron, i.e. includes the kinetic energy. The 
Monte-Carlo data, if plotted, would be indistinguishable on 
the scale used here from the curves shown. For reference, 
we also show the total energy per electron (dotted line) ac- 
cording to an alternative parameterisation, the well-known 
Wigner interpolation formula, Ref. [K2|. 



C. The Local Density Approximation 

As a practical approximate expression for E^^ ['^•1 , Kohn 
and Sham p7[ | suggested what is known in the context of 
DFT as the local density approximation, or LDA: 



£;xc[n(r)] 



dv n(r) exc(ri(r)) , 



(22) 



where exc('^) is the exchange-correlation energy per elec- 
tron in a uniform electron gas of density n. This quan- 
tity is known exactly in the limit of high density, and 
can be computed accurately at densities of interest, us- 
ing Monte Carlo techniques (i.e. there are no free pa- 
rameters). In practice one usually employs paramet- 
ric formulas, which are fitted to the data and are ac- 
curate to within 1-2%. As an example, we quote that 
given by Gunnarson and Lundqvist (Ref. [^l[): exc("-) — 
-0.458/rs - 0.0666G'(rs/11.4) Hartrees, where is in 
units of the Bohr radius, and G{x) = i{(l + a;'^) log(l-f' 
x-i) -x^ + ^x- i}, see Fig. |. 

Note that the only difference between the resulting 
computational scheme and a naive mean-field approach 
is the addition of the potential 



i'xc(r) 



d (nebcin)) 



dn 



(23) 



to the electrostatic potential at the appropriate step in 
the self-consistency loop p^ ] . The corresponding expres- 
sion for the ground-state energy is: 

N 

Eo = ^ei- Ees[n{r)] + dr n(r) (exc(n(r)) - v^c{n{r))) 

i=l ■' 

(24) 

where the first term is the noninteracting energy, the sec- 
ond term subtracts half of the double counting of the elec- 
trostatic energy as in the Hartree scheme, and the last 
term is a similar subtraction for the exchange-correlation 
energy. 

The LDA has been shown to give very good results for 
many atomic, molecular and crystalline interacting elec- 
tron systems, even though in these systems the density 
of electrons is not slowly varying. As an example, we 
show in Fig. ^ the solution for an atom of Argon. One 
can see that the shell structure, which is absent in the 
Thomas-Fermi model, is described here in detail. The 
calculated ground-state energy is —525.9 Hartrees, com- 
pared to —652.7 Hartrees in the Thomas-Fermi model 
(one Hartree = me^/h^ = 27.2eV). The experimental 
value is —527.6 Hartrees — the LDA result is accurate 
to within less than half a percent, compared to the 20% 
inaccuracy of the Thomas-Fermi model. Both the de- 
scription of the shell structure and much of the improve- 
ment in the energy estimate are due to the introduction 
by Kohn and Sham of an exact method for evaluating 
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Ar Z=18 




dA 



(25) 



FIG. 4. The density n(r) for an Argon atom, multiplied 
by for convenience of presentation (inverse Bohr radii), as 
a function of the distance from the origin (in Bohr radii). The 
full line is the DFT result using the LDA, and shows the shell 
structure; the dashed line is the result of the Thomas-Fermi 
model for this atom. 



the highly nonlocal kinetic energy functional. However, 
the accuracy of the exchange-correlation term is also of 
central importance, and will be discussed next. 



IV. ACCURACY OF THE 
EXCHANGE-CORRELATION ENERGY 

In this section, we discuss a few aspects of i?xc [n] which 
are somewhat more advanced (and may be skipped in a 
first reading). In the first subsection, a formally exact 
expression for the exchange-correlation energy is derived, 
and a sum-rule which applies to it is obtained. The next 
subsection provides references to some more recent and 
accurate approximations, which were devised with this 
sum-rule in mind. Finally, a general discussion of the 
level of accuracy achieved by these approximations for 
different applications is given. 



A. The exchange— correlation hole 

A deeper understanding of the exchange-correlation 
energy can be achieved by considering a continuous tran- 
sition between the interacting and noninteracting sys- 
tems which appear in the definition of i?xc[?i], rather 
than using a simple subtraction as in Eq. disl). To do 
this, we reduce the strength of the Coulomb interaction, 
e^, and contemplate a more general interaction potential, 
u{r) = Ke^ /r with < A < 1. In other words, we de- 
fine the Hamiltonian as Hmb — T + V + AU, with the 
noninteracting system corresponding to A = 0, and the 
intcracting-electron system to A = 1. For a given density 
distribution n(r) one can consider i^HKi^^] for intermedi- 
ate values of A, which leads to the exact expression [k4[ : 



The general rules for Legendre transforms give the deriva- 
tive (9FhkM/9A as equal to dil/dA, which in turn is 

equal to (U), the expectation value of the interaction 



i?xcW = /odA 
correlation energy 



energy U. Note that the derivative of fl is taken at a 
constant potential, but in order to reproduce the given 
density distribution n{r), this potential must depend on 
A (it is usually denoted by t;A(r), but we need it here 
only at the two extremes, for which we already have ap- 
propriate notations: wi(r) — v{r) and uo(r) — Wcff(r)). 
Comparing Eq. ( p5| ) with Eq. (p^, we find that 

UJ — £'cs[^], i.e. the exchange- 
is the difference between the A- 
averaged expectation value of the interaction energy and 
the electrostatic approximation to it. Both of these are 
integrals over the product of the Coulomb interaction and 
the density of pairs, which is ^ (/5(r) (jo(r') — 5{r — r'))) 
and in(r)n(r') for the exact and approximate expres- 
sions, respectively (the factor of ^ corrects for the dou- 
ble counting of each pair, and the S function removes 
the interaction of each electron with itself). We are thus 
motivated to define a pair correlation function, g{r,r'), 
equal to the difference between the two pair densities 

ff(r,r') = (/5(r)(p(r')-<S(r-r')))„(r),A" ^W«(r')- 

(26) 

It is of relevance here that g{r, r') at typical electron den- 
sities is relatively featureless — it is only at very low den- 
sities {ts ~ lOOao) that electrons tend to develop strong 
correlations, and may form a Wigner crystal. 

Using the fact that the number of particles, M = 
J dr' /3(r'), does not fiuctuate at very small temperatures, 
and is equal to its expectation value, N = J dr' n(r'), it is 
easy to show that / dr'g(r,r') — —n{r). It is convenient 
to define a normalized version of this correlation func- 
tion pxc(r,r'; [n]. A) = g{r,r')/n{r), which is the density 
of the so-called exchange-correlation hole, pxc describes 
the region in r'-space from which an electron is "miss- 
ing" if it is known to be at the point r, and the fact that 
its integral over r' is equal to — 1 corresponds to the fact 
that there is exactly one "missing" electron. In terms of 
the A-integrated value of this quantity, 



w(r, r' 



')= /'dApxc(r,r';[n],A)= dA ^(r, r')/n(r) 
Jo Jo 



(27) 



the exchange-correlation energy is given formally as 
f drdr' , 



(28) 
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In other words, the exact exchange-correlation energy 
may be written as in the LDA, Eq. (|2|), provided 
that the corresponding energy per particle, Exc is in- 
terpreted not as a local quantity but is evaluated ac- 
cording to the density of the exchange-correlation hole, 
exe = (eV2)/dr'«;(r,r')/|r-r'|. 

The normalization prop- 

erty of the exchange-correlation hole, pxc, carries over 
to its A-averaged counterpart: 



J dr' 'ii;(r,r') = -1 



(29) 



This "sum-rule" , Eq. ( P9| ) , has been used Q as the ba- 
sis for an "explanation" of the relatively high accuracy 
achieved by the LDA: the reference system here (the ho- 
mogeneous electron gas) has properties which are also 
exact for the inhomogeneous system. More importantly, 
it restricts and guides the search for more accurate prac- 
tical approximations: expressions which break the sum- 
rule can not be expected to work well (see below). 



B. Refinements of £xc 

The fact that the LDA achieves a high relative accu- 
racy, i.e. can predict the ground-state energy of various 
systems to within less than a percent, does not mean 
that the absolute accuracy is sufficient. In particular, 
typical applications in chemistry require that the en- 
ergy of a molecule be known to within a small fraction 
of an electron-volt. Indeed, improving upon the accu- 
racy of the LDA is a goal which has been persistently 
pursued. One improvement which is very often imple- 
mented is the local spin-density (LSD) approximation 
p6| , which is motivated in part by the fact that the 
exchange-correlation hole is very different for electrons 
with parallel and with antiparallel spins. In this scheme, 
separate densities of spin-up and spin-down electrons are 
used as a pair of functional variables: ni{r) and n^lr), 
and the Hamiltonian contains separate potentials for spin 
up and spin down electrons — a Zeeman-energy mag- 
netic field term is introduced. The exchange-correlation 
energy per particle is then taken from the results for 
a homogeneous spin-polarized electron gas, exdni^ni). 
The spin dependence allows Hund's rule to be discussed 
within DFT. 

The next degree of sophistication is to allow e^c to 
depend not only on the local densities but also on the 
rate-of-changc of the densities, i.e. to add gradient cor- 
rections. Unfortunately, it was found that such correc- 
tions do not necessarily improve the accuracy obtained. 
In fact, introducing gradient corrections in a straightfor- 
ward and systematic manner, by expanding around the 
uniform electron gas, breaks the sum rule of Eq. ( |29| ) 
and is less accurate |^]. This situation led to the de- 
velopment of various generalized gradient approximations 



(GGAs) |^,|8[, in which the spatial variations of n(r) en- 
ter in a manner which conforms with the sum rule, and 
which have succeeded in reducing the errors of the LDA 
by a factor which is typically about 4. 

Further improvements in practical expressions for 
Exc[n\ are actively being pursued j3|]. One direction 
which may perhaps achieve the accuracy needed for ap- 
plications in chemistry [^, is to use the fact that the 
exact form of the exchange-correlation hole can be cal- 
culated for A = relatively easily, directly from the non- 
interacting Kohn-Sham system. There is thus no need 
to use an approximation such as the LDA or the GGA 
for the low-A portion of the integral in Eq. (^). Ulti- 
mately, one hopes that a systematic method of improving 
the approximation would be found, although so far this 
has been an elusive goal. 



C. Successes and failures 

Over the years, many different types of applications 
of DFT have been developed. This variety evolved be- 
cause knowledge of the electronic ground-state energy as 
a function of the position of the atomic nuclei determines 
molecular and crystal structure, and gives the forces act- 
ing on the atomic nuclei when they are not at their equi- 
librium positions. At present, DFT is being used rou- 
tinely to solve problems in atomic and molecular physics. 



such as the calculation of ionization potentials |41| and 
vibration spectra, the study of chemical reactions, the 
structure of bio-molecules , and the nature of active 
sites in catalysts [ p3| , as well as problems in condensed 
matter physics, such as lattice structure s , phase tran- 
sitions in solids [^, and liquid metals |46[. Furthermore 
these methods have made possible the development of ac- 
curate molecular dynamics schemes in which the forces 
are evaluated quantum mechanically "on the fly" ji^ j. 

It is important to stress that all practical applications 
of DFT rest on essentially uncontrolled approximations, 
such as the LDA discussed above. Thus the validity of 
the method is in practice established by its ability to 
reproduce experimental results. A discussion of the ac- 
curacy achieved by DFT, compared to other alternative 
approaches, necessarily depends very much on the spe- 
cific applications one has in mind, as detailed below. 

For atoms and small molecules, the simplest version 
of the LDA already provides a very useful qualitative 
and semi-quantitative picture. It is of course a dra- 
matic improvement over the Thomas-Fermi model. It 
even improves on the more labor-intensive Hartree-Fock 
method in many cases, especially when one is calculat- 
ing the strength of molecular bonds, which are substan- 
tially overestimated in Hartree-Fock calculations. This 
can only be considered as a surprising success, keeping in 
mind that an isolated atom or molecule is as inhomoge- 
neous an electronic system as possible, and therefore the 
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last place where one might expect a local approximation 
to work. In other words, electronic correlations in such 
systems are in a sense weak, and are on average similar 
to those of a uniform electron gas (see the discussion of 
the sum-rule, Eq. (p9|)). However, the many-body quan- 
tum states of such relatively small systems can be solved 
for extremely accurately using well-known techniques of 
quantum chemistry, specifically the configuration inter- 
action (CI) method Q. Furthermore, these techniques 
use controlled approximations, so that the accuracy can 
be improved indefinitely, given a powerful enough com- 
puter, and indeed impressive agreement with experiment 
is routinely achieved. For this reason, most quantum 
chemists did not embrace DFT at an early stage. 

It is in studies of larger molecules that DFT be- 
comes an indispensable tool The computational ef- 
fort required in the conventional quantum chemistry ap- 
proaches grows exponentially with the number of elec- 
trons involved, whereas in DFT it grows roughly as the 
third power of this number. In practice, this means that 
DFT can be applied to molecules with hundreds of atoms, 
whereas using CI, one is limited to systems with only a 
few atoms. Simply solving the noninteracting problem 
for a complicated molecule may also be prohibitive, and 
various methods are used in order to reduce the problem 
to a computationally manageable task. Of these, we men- 
tion the well-known pseudopotential method , which 
allows one to avoid recalculating the wavefunctions of 
the inert core electrons over and over again, and the re- 
cent attempts to develop "order N" methods which 
make use of the fact that the behavior of the densities 
at each point is determined primarily by the atoms in its 
immediate vicinity, rather than by the whole molecule. 
It is for this problem that more and more accurate den- 
sity functionals are most obviously needed. To illustrate 
this, we quote one sentence from Ref. js^ j: "Accurate at- 
omization energies are found [using the GGA] for seven 
hydrocarbon molecules, with a rms error per bond of 0.1 
eV, compared with 0.7 eV for the LSD approximation 
and 2.4 eV for the Hartree-Fock approximation." 

The remarkable usefulness of DFT for solid-state 
physics was apparent from the outset. For example, the 
lattice constants of simple crystals are obtained with an 
accuracy of about 1% already in the LDA |^^. In such 
applications, the electronic structure of a single unit cell 
with periodic boundary conditions is studied; more ambi- 
tious applications are also common, e.g. a supercell con- 
taining many unit cells with a single impurity or defect 
pTf . Admittedly, this method is inappropriate for treat- 
ing some more complicated situations, such as antiferro- 
magnets or systems with strong electronic correlations. 
In other cases, such as for the work-function of metals, 
local approximations such as the LDA obviously miss an 
important part of the physics: for a point r a short dis- 
tance away from the surface of a metal, the exchange- 
correlation hole pxc(r, r') is concentrated at points r' in- 



side or very near the surface of the metal; this results in 
image forces, i.e. a l/r behavior of Vxc (where r is the 
distance from the surface) which is nonlocal. However, 
this deficiency can be corrected for "by hand" , yielding 
satisfactory results . 

In general it is useful to note that, in contrast to ap- 
proximations using free parameters which are empirically 
optimized to fit a certain set of data and may thus be used 
reliably for interpolation, the LDA and the GGA have 
proved to exhibit a consistent degree of accuracy or inac- 
curacy for a wide variety of problems — when applied to 
a new problem, the results can thus be interpreted with 
some confidence. One should of course also be aware of 
the cases for which these approximations are known to 
fail, such as the image forces mentioned above, and van 
der Waals forces |53|| , which are important e.g. for bi- 
ological molecules. Both of these are manifestations of 
the significance of nonlocal correlations — a nonlocality 
which is by definition absent from the LDA and its im- 
mediate extensions. These examples of practical failure, 
together with the unattractiveness of uncontrolled ap- 
proximations, spur research towards new and more exact 
exchange-correlation energy functionals. 

Our discussion would not be complete without men- 
tioning the existence of many other uses of density- 
functional methods, for electronic systems and for other 
physical systems. The former include time-dependent 
DFT, which relates interacting and noninteracting elec- 
tronic systems moving in time-dependent potentials, and 
relativistic DFT, which uses the Dirac equation rather 
than the Schrodinger equation to calculate the Kohn- 
Sham states (these are reviewed in Ref. [^). The latter 
include applications in nuclear physics, in which the den- 
sities of protons and neutrons and the resulting energies 
are studied [Q, and in the theory of liquids, as already 
discussed in Sec. II. 

V. SUMMARY 

In describing density functional theory (DFT) and the 
approximations typically implied by its use, it is nec- 
essary to follow two steps, as was done in considerable 
detail in Sees. II and HI above. The first step is to in- 
troduce the Hohenberg-Kohn energy, FHK[?i(r)]. It is 
equal to the internal energy, i.e. the difference between 
the ground-state energy Eq and the potential energy, or 
the sum of the kinetic energy and the interaction energy 
of the electrons. i^HKl'T-] is a universal functional of the 
density distribution — it applies to atoms, molecules, 
crystals, and all other electronic systems. The existence 
of i^HK["-] arises from the fact that each system has not 
only a unique external potential, v{r), as in traditional 
many-body theory, but also a unique density distribu- 
tion, n{r). Within DFT, the different systems are la- 
beled by their different electronic densities, n{r), and the 
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potential v{r) is considered as secondary to, and depen- 
dent on, the primary n(r). In fact, the potential is given 
by the functional derivative v{r) — —6FHK/Sn{r). This 
property of FnKi'T-] can be deduced in two equivalent 
ways: (a) it is the Euler equation for the Hohenberg- 
Kohn theorem, which is a variational principle stating 
that the total energy, Fhk["] + / dr n{r)v{r), is min- 
imized for a given v{r) by the corresponding ground- 
state density n(r), and has the ground-state energy Eo 
as its minimum value; (b) it arises as the conjugate of 
the well-known relationship n(r) = SEo/6v{r), when the 
definition FHRi^T-] = Eq — J dr n{r)v{r) is viewed as a 
(fimctional) Legendre transform. 

The functional FHK[?^(r)] discussed above is not known 
explicitly in terms of its variable, n{r). The second cru- 
cial step is to introduce a practical approximation for 
this energy functional. It may be written as a sum of 
three terms: (i) the kinetic energy term, which is by 
definition equal to the value FHK["-(r)] would have for 
noninteracting electrons; (ii) the Hartree term, which is 
simply a double integral over n(r)n(r'); and (iii) the re- 
mainder or exchange-correlation term, which in princi- 
ple contains all of the complicated interaction physics 
ignored by the first two terms, and must be approxi- 
mated in practice. By taking a functional derivative, 
one finds that the external potential v{r) is given by a 
corresponding sum of three terms: (i) the so-called ef- 
fective potential, Wcff(r), which is the potential which a 
system of noninteracting electrons must have in order 
to reproduce the 7i(r) of the interacting system, (ii) the 
electrostatic potential, and (iii) an exchange-correlation 
potential. The simplest widely-used expression for the 
exchange-correlation term is the local density approx- 
imation (LDA), which takes the exchange-correlation 
energy-density at each point in the system to be equal 
to its known value for a uniform interacting electron gas 
of the same density, and results in a parameter-free ap- 
proximate description of all electronic systems. 

These ideas are implemented by the Kohn-Sham set of 
equations, which consists of a noninteracting Schrodinger 
equation involving the effective potential Woff(r), the 
abovementioned relationship between Wcff(r) and the 
given external potential u(r), and expressions for the den- 
sity distribution n{r) and the interacting ground-state 
energy Eq in terms of the properties of the noninteracting 
single-particle solutions and of the approximate expres- 
sion for the exchange-correlation energy which is in use. 
Using modern computers, the Kohn-Sham equations can 
be solved even for systems containing dozens of atoms, 
and the results for the ground-state energy are typi- 
cally accurate to within a small fraction of a percent. In 
contrast to other (much more computationally demand- 
ing) methods of calculating the ground-state energy, the 
LDA is an uncontrolled approximation, and so there is no 
straightforward path to desired further improvements in 
the accuracy. Nevertheless, remarkable progress in this 



direction has been achieved over the years, most notably 
with the introduction of the generalized gradient approx- 
imation (GGA). 

Whereas the outline just given could apply (with mi- 
nor modifications) to other introductions to DFT, the 
present discussion was based on an analogy with ther- 
modynamics. It is well known that in treating situations 
where the the number of particles N is constrained, it 
is preferable to use it as a free variable, rather than the 
chemical potential /i. Similarly, one may think of the 
Coulomb interaction as imposing strong constraints on 
the density distribution n(r) required to achieve low- 
energy structures in inhomogeneous electronic systems; 
it is thus preferable to use it instead of the potential v{r) 
as a free variable. 

Three of the advantages of the present approach, 
as compared, e.g., with introducing DFT using Levy's 
constrained-search method [Q, are: (a) the density dis- 
tribution n(r) appears here as a natural variable — it 
is conjugate to v{r) through a Legendre transform — 
whereas in the conventional description of DFT the very 
existence of the functional Fhk [n] appears to be surpris- 
ing and requires some digestion; (b) some of the mathe- 
matical difficulties encountered in the ground-state the- 
ory are not present in the theory of finite temperature 
ensembles; and (c) using the standard properties of Leg- 
endre transforms, one immediately obtains the physical 
expression for the exchange-correlation energy in terms 
of the density of the exchange-correlation hole, Eq. ([2^), 
an expression which serves as the basis for a discussion 
of the weaknesses and strengths of the approximations 
employed in practice. We hope that the availability of 
this type of introduction will help increase the awareness 
and understanding of DFT amongst potential users, and 
especially amongst the general audience of physicists and 
scientists. 
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